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Abstract 

In this paper we consider generalized nonlinear Schrodinger equations with ex- 
ternal potentials, we compute the forth and the sixed order Fermi Golden Rules 
(FGR), conjectured in GS2l IGS3) . which is used in a study of the asymptotic 
dynamics of trapped solitons. 

1 Introduction 

It is well known that the eigenvalues embedded in a continuous spectrum of 
a self-adjoint operator, L, are unstable under generic perturbations. Intu- 
itively, one can think of time periodic solutions ('bound states') of the equa- 
tion, iii = Lu, corresponding to such eigenvalues, leaking their energy to the 
continuous spectrum solutions ('scattering states') which in turn is radiated to 
infinity. The coupling between time-periodic and continuous spectrum solutions 
was computed in the second order of perturbation theory by P. A. M. Dirac in 
1927 (see |Dij . resulting in an elegant expression. If this expression is nonzero, 
then the decay mentioned above takes place. E. Fermi called this condition 
"Golden Rule No. 2" . In physics literature it is known as the Fermi Golden 
Rule (FGR). The FGR, computed first for atomic and molecular states, appears 
in many other areas of physics e.g. in non-equilibrium statistical mechanics. Of 
course, due to energy conservation the isolated eigenvalues are stable under 
reasonable time independent perturbations. 

The Fermi Golden Rule was introduced to nonlinear Hamiltonian PDE in S 
where it was used to prove that time-periodic solutions for linear (and nonlinear) 
wave equations are unstable under generic nonlinear perturbations. It was used 
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extensively in investigating long-time dynamics of solitary waves in the nonlinear 
Schrodinger |BF2I IBlISul ISWH ISW21 IBW5I IT7T1 IT72I IT75| and wave fSW4| 
equations. 

However, in nonlinear problems with a small parameter, h, of the quasi- 
classical nature, the second order expression obtained by Dirac vanishes and, in 
principle, one can hope to get non-vanishing expression only in the 2A^th order 
where N = O(j^). (If one thinks about the FGR as a condition that an 'atom is 
ionized by a single photon', in this case because the ionization energy is much 
larger than the photon frequency, this 'ionization' is due to a multi-photon 
process.) This type of situation occurs in particular in the problem of the 
relaxation of a solution to the ground state in generalized nonlinear Schrodinger 
equations with external slowly varying potentials (Gross-Pitaevskii equations) 
described below. 

In this paper we find, for the first time, the 'nonlinear' forth and sixth order 
FGRs for the relaxation problem mentioned above. As was conjectured in GS2 
GS3] , they are of the same form as the second order FGR and as, presumably, 
the higher order FGRs. To obtain the forth and sixth order FGRs we have 
to perform additional normal-form-type transformations in the perturbation 
theory we developed in |GS2llU5S] . 

Problem. Now we describe the problem in more details. We consider the 
generalized nonlinear Schrodinger equation (NLS) in dimension d ^ 2 with an 
external potential Vh '■ K d — > K, 

i^ = -A^ + V h ^-fM 2 )^. (1) 

Here h > is a small parameter giving the length scale of the external potential 
in relation to the length scale of the Vh = solitons (see below) , A is the Laplace 
operator and /(s) is a nonlinearity to be specified later. We normalize /(0) = 0. 
Such equations arise in the theory of Bose-Einstein condensation , nonlinear 
optics, theory of water waves 2 and in other areas. 

To fix ideas we assume the potentials to be of the form Vh(x) :— V(hx) 
with V smooth and decaying at 00. Thus for h — 0, Equation (QJ becomes the 
standard generalized nonlinear Schrodinger equation (gNLS) 

^ = -A^ + / ^-/(b/f)V, ( 2 ) 

where fi := V^(0). For a certain class of nonlinearities, f{\ip\ 2 ), ( see Section |3J), 
there is an interval Tq C K such that for any A S Tq Equation ( [2J has solutions 
of the form e'^^Vo ( x ) where 4>^ e W 2 (R ra ) and ^ > 0. Such solutions (in 
general without the restriction c/>g > 0) are called the solitary waves or solitons 
or, to emphasize the property > 0, the ground states. For brevity we will use 
the term soliton applying it also to the function </>q without the phase factor 

e i(A-jLi)t_ 

1 In this case Equation (0 is called the Gross-Pitaevskii equation. 

2 In the last two areas Vh arises if one takes into account impurities and/or variations in 
geometry of the medium and is, in general, timc-dcpcndcnt. 
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Equation (HJ is translationally and gauge invariant. Hence if e^ A ^^(x) 
is a solution for Equation (|3J), then so is 

e^-^e^^ix + a), for any a S R", and a G [0, 2tt). 

This situation changes dramatically when the potential Vh is turned on. In 
general, as was shown in |F Wl lObiTl lABCj out of the (n + 2)-parameter family 
e j(A-p)t e ia0A^ x _j_ a ^ Qm y a discrete get f two-parameter families of solutions 
to Equation ([TJ bifurcate: e 4At e M A (x), a S [0, 2tt) and A e I for some 1 C 
2"o> with <f) X = <fi^ E H.2 (K") and A > 0. Each such family centers near a 
different critical point of the potential Vh(x). It was shown in |Oh2| that the 
solutions corresponding to minima of Vh{x) are orbitally (Lyapunov) stable and 
to maxima, orbitally unstable. We call the solitary wave solutions described 
above which correspond to the minima of Vh{x) trapped solitons or just solitons 
of Equation ( ^) omitting the last qualifier if it is clear which equation we are 
dealing with. 

A basic question about soliton solutions is whether they are asymptotically 
stable, i.e. whether for initial condition of ( ^| sufficiently close to a trapped 
soliton {e n <p x (xj) the solution converges (relaxes) in a local norm, up to a phase 
factor, to another trapped soliton, 

V>(z,i)-e i7(t) Aoo (:r) ^0. 

We observe that ( ^ is a Hamiltonian system with conserved energy (see 
Section |2J and, though orbital (Lyapunov) stability is expected, the asymptotic 
stability is a subtle matter. To have asymptotic stability the system should be 
able to dispose of excess of its energy, in our case, by radiating it to infinity. 
The infinite dimensionality of a Hamiltonian system in question plays a crucial 
role here. This phenomenon as well as a general class of classical and quantum 
relaxation problems was pointed out by J. Frbhlich and T. Spencer [Private 
Communication] . 

Another important property is their effective dynamics. Namely, one would 
like to show that if an initial condition is close (in the weighted norm :— 
||(1 + |a;| 2 )^M| |^i for sufficiently large v) to the soliton e ll0 4> x ° , with -f g K and 
A G T (X as above), then the solution, x/)(t), of Equation (^| can be written as 

ijj{x,t) = e ll(t) [e lp(t) - x (j) X{t) {x-a{t))+R{x,t)]. / (3) 

where |ji?(i)||_„,i — > 0, A(i) —> Aoo for some Aoo as t — * oo and the soliton 
center a(i) and momentum p(t) evolve according to certain effective equations 
of motion. 

Results. As in |GS21 lUS^j we assume that either d = 1 and the potential 
is even or d > 2 and the potential is spherically symmetric and the initial 
condition symmetric with respect to permutations of the coordinates. In this 
case the soliton relaxes to the ground state along the radial direction. This 
limits the number of technical difficulties we have to deal with. We expect that 
our techniques extend to the general case when the soliton spirals toward its 
equilibrium. 
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It is shown in |GS21 lUSlj] . under certain conditions, that the ground state 
is asymptotically stable and the effective equations for the parameters a and p 
are close to Newton's equations. It was conjectured in |GS2I IGS3| that one of 
the conditions (which plays a key role) is equivalent to an explicit 2iVth order 
FGR. This conjectured is tru e for N = 1 due to the works [BPll IBuSul ISWT1 
[?TW^I?W^lT^lTY2llT73| . In this paper we prove this conjecture for N = 2 
and 3. In other words we establish the FGR for nonlinear problem in the 4th 
and 6th orders. 

Moreover, we find more precise effective equations for a and p. We show 
that there exist a function z(t) and a constant Zn+i,n € C satisfying z(t) = 
p{t) + ia(t) + 0(\p\ 2 + \a\ 2 ) and 

lf t \z\ 2 =ReZ N+l , N \z\ 2N+2 + 0(\z\ 2N +*) 

with N — 0(1 /h) being an integer depends on A and h. For N = 2, 3 we find 
the explicit form of ReZff+i,N which in particular shows that it is always non- 
positive and negative generically. The explicit form is given by 2iVth order FGR 
mentioned above. 

In the present paper by using normal forms we simplify the expressions 
for z and the other parameters, thus make the proof more transparent than 
|HS2llHS3l . 

Previous results. We refer to |GS1| for a detailed review of the related 
litera t ure. Here we only mention results of [Cu| IEFT1 lBF2l IBuSul ISWT1 ISWS1 
ISW51 ITYT1 ITY21 ITY5] which deal with a similar problem. Like our work, 
|Wl1IW2lRWlllTYl1lTT2llT73] study the ground state of the NLS with a 
potential. However, these papers deal with the near-linear regime in which the 
nonlinear ground state is a bifurcation of the ground state for the corresponding 
Schrodinger operator —A + V(x). The present paper covers highly nonlinear 
regime in which the ground state is produced by the nonlinearity (our analysis 
simplifies considerably in the near- linear case). 

Papers [Cu| |BP"T1 IEF21 IBuSuj consider the NLS without a potential so 
the corresponding solitons, which were described above, are affected only by a 
perturbation of the initial conditions which disperses with time leaving them 
free. In our case they, in addition, are under the influence of the potential and 
they relax to an equilibrium state near a local minimum of the potential. 

As was mentioned above, the relaxation of the solution to the ground state 
(the asymptotic stability of the ground state) was shown in |GS2I [GS3j under 
some assumptions on the potential V , nonlinearity / and the linearized operator 
L(X), which in particular include these related to the FGR which was explained 
above. 

The paper is organized in the following way. In Section \T\we introduce the 
concept of Fermi Golden Rules and outline its applications and importance in 
proving the asymptotic stability of trapped solitons of nonlinear Schrodinger 
equations. In Section [2] we show the conservation laws and the local well- 
posedness of (0. In Section |3|we formulate the conditions on the nonlinearity 
/ and the potential V so that there exists a soliton manifold to ( and it is 
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stable. In Section |||we linearize the solution to (QJ around the soliton get a 
linear operator. Moreover we analyze its spectrum and put some assumptions. 
In Section [SJwe state the main theorem. In Section Elwe separate the 'useless' 
parts from the equations of parameters, to prepare for the computation of forth 
(N = 2) and sixth (N = 3) order Fermi Golden Rules which will be proved in 
Section 0and |S] respectively. 

Notation. As customary we often denote derivatives by subindices as in 
(f)\ = -§^4> x for <j) X = 4> x (x). However, the subindex h signifies always the de- 
pendence on the parameter h and not the derivatives in h. The Sobolev and L 2 
spaces are denoted by H k and C 2 respectively. 

Acknowledgment. This paper is part of the author's Ph.D thesis require- 
ment. The author wishes to thank his advisor, I.M.Sigal, for introducing the 
problem and useful discussions. 

We are grateful to J. Colliander, S. Cuccagna, S. Dejak, J. Frohlich, Z. Hu, 
W. Schlag, A. Soffcr, G. Zhang, V. Vougalter and, especially, V.S. Buslaev for 
fruitful discussions. Upon the completion of this paper we learn from T.P.Tsai 
that he also obtained the explicit form of Fermi Golden rule for N = 2 [Private 
Communications] . 

2 Hamiltonian Structure and GWP 

Equation (QJ is a Hamiltonian system on Sobolev space W 1 (M, C) viewed as a 
real space H 1 (R, R) ® with the inner product (tp, <f>) — Re Jr?/^ and 

with the simpletic form Lu(ip,(ft) — Im J R ip<p. The Hamiltonian functional is: 

m) ■■= f[\(m 2 + v h \M 2 )~F{\i>\ 2 % 

where F(u) :— \ f™ f(£)d£. 

Equation ( ^| has the time-translational and gauge symmetries which imply 
the following conservation laws: for any t > 0, we have 

(CE) conservation of energy: H(ip{t)) = H(ip(0)); 

(CP) conservation of the number of particles: N(ifj(t)) = N(ip(0)), where N(tp) := 

To address the global well-posedness of ( 0) we need the following condition on 
the nonlinearity /. Below, s + — s if s > and = if s < 0. 

(fA) The nonlinearity / satisfies the estimate |/ (£)| < c(l + |<^| Q_1 ) for some 
a e P. (3^) and 1/(01 < "(1 + leP) for some /3 e [0, |). 

The following result can be found in |Caz| . 
Theorem Assume that the nonlinearity f satisfies the condition (fA), and that 
the potential V is bounded. Then Equation is globally well posed in Ji , i.e. 
the Cauchy problem for Equation fQ|) with a datum ip(0) £ Ji 1 has a unique 
solution i/j(t) in the space TL l and this solution depends continuously on ip(0). 
Moreover tj){t) satisfies the conservation laws (CE) and (CP). 
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3 Existence and Orbital Stability of Solitons 

In this section we review the question of existence of the solitons (ground states) 
for Equation ( ^) . Assume the nonlinearity / : M — > R is smooth and satisfies 

f(s) 

(fB) There exists an interval^ G ^ + s -t- for any A G Tq, — oo < lim s _, +00 — -y- < 
and | J* f(s)ds > A for some constant £, for d > 2; and 

C/(0,A) :=-A0 2 + / 

has a smallest positive root </>o(^) such that [7^(00 (A), A) ^ 0, for d = 1. 

It is shown in |BLI IStr| that under Condition (fB) there exists a spherical 
symmetric positive solution tfi x to the equation 

-A0 A + \<j) X - f(^ x ) 2 )0 x = 0. (4) 

Remark 1. Existence of soliton functions <f> x for d — 2 is proved in tStrf under 
different conditions on f. 

When the potential V is present, then some of the solitons above bifurcate 
into solitons for Equation (2J. Namely, let, in addition, / satisfy the condition 
1/ (C)l — c (l + IC| p )i for some p < oo, and V satisfy the condition 

(VA) V is smooth and is a non-degenerate local minimum of V. 

Then, similarly as in jFWI IOhl| one can show that if h is sufficiently small, 
then for any A € Zqv, where 

1 QV := {A|A > - inf {V(x)}} f] {A|A + V(0) € Z }, 

there exists a unique soliton </> A = (\>\ (i. e. (f> x G Ti.2(R) and <fi x > 0) satisfying 
the equation 

-A0 A + (A + v h )4> x - f{^ x f)4> x = o 

and the estimate <f> x — <fi X+V ^ ' + 0(h 3 ^ 2 ) where <f> x is the soliton of Equation 
(0. 

Let 8{X) := \\\4> x \\\- It is shown in |QSS1| that the soliton 4> x is a minimizer 
of the energy functional H(tp) for a fixed number of particles N(^) = constant 
if and only if S (A) > 0. Moreover, it is shown in (We2l IGSS1| that under 
the latter condition the solitary wave tfi x e lXt is orbitally stable. Under more 
restrictive conditions (see |CjSS1| ') on / one can show that the open set 

1 := {A S lov : S'(X) > 0} (5) 

is non-empty. Instead of formulating these conditions we assume in what follows 
that the open set X is non-empty and A G T. 

Using the equation for 4> x one can show that if the potential V is re- 
dially symmetric then there exist constants c, 5 > such that |0 A (x)| < 
ce -5 M and |^^ A | < ce - * 5 ' 1 ', and similarly for the derivatives of 4> x and -j^4> X - 
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4 Linearized Equation and Resonances 

We rewrite Equation ([IJ as = G(ip) where the nonlinear map G(tp) is defined 
by G(ip) = — i{— A + A + Vh)ip + if(\ip\ 2 )ip. Then the linearization of Equation 
([TJ can be written as ^ = dG(4> x )x where dG((f> x ) is the Frechet derivative of 
G(ip) at <f>. It is computed to be 

dG(cp x ) X = -t(-A + A + + */(W> A ) 2 )x + 2 J /'((0 A ) 2 )(0 A ) 2 i? ex . (6) 

This is a real linear but not complex linear operator. To convert it to a linear 
operator we pass from complex functions to real vector-functions x < — * X — 

^ ) , where xi = R e X an d Xi = Imx- Then dG((p x )x < — ► L(X)x where 
X2 J 

the operator L(X) is given by 

m - ( \ m fr m ) • P) 

with L_(A) :--A + 14 + A-/((0 A ) 2 ),and L+(A) := — A + V/j + A — f((4> x ) 2 ) — 
2/'(O A ) 2 )(0 A ) 2 . The operator L(A) is extended to the complex space H 2 (K, C)© 
7i 2 (R, C). If the potential Vh in Equation decays at oo, then by a general 
result 

er ess (£(A)) = (-ioo, —iX] n [iA, zoo). 
The eigenfunctions of £(A) are described in the following theorem (cf GS1 . 

ESU). 

Theorem 4.1. Let V satisfy Condition (VA) and \h\ be sufficiently small. Then 
the operator L(X) has at least 2d + 2 eigenvectors and associated eigenvectors 
with eigenvalues near zero: two-dimensional space with the eigenvalue and a 
2d- dimensional space with non-zero imaginary eigenvalues ±iej(X), 

ej (X) := hy/2e~j + o(h), 

where ej are eigenvalues of the Hessian matrix ofV at value x = 0, V (0). The 
corresponding eigenfunctions ( J are re ^ ec ^ by complex conjugation and 

satisfy 

£j = \/2d Xk (j) X + o(h) and r/j = —hy/lijXj4> X + o(h), 
and £i and rjj are real. 

Remark 2. The zero eigenvector ^ ^ A and the associated zero eigenvector 
(b x \ 

q J are related to the gauge symmetry t) ~ * e la ip(x, t) of the original 
equation and the 2d eigenvectors ( j!^. 1 with 0(h) eigenvalues originate 
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(d d x \ 
Q Xk 1 , A; = 1, 2, • • -, d, and the associated zero 

eigenvectors ^ ^ ^ A , k = 1, 2, • • •, d, of the V — equation due to the 

translational symmetry and to the boost transformation ip(x,t) — > e lb ' x ip(x,t) 
(coming from the Galilean symmetry), respectively. 

For d > 2 we will be interested in permutationally symmetric functions, 
g G C 2 (R d ), characterized as 

g{x) = g(crx) for any a E S d 

with Sd being the group of permutation of d indices and a(x\,X2,- ■ -,Xd) ■= 
{x a (i),x a{ 2), ■ ■ -,a;<T(d))- 

Remark 3. For any function of the form e tp ' x <p(\x — a\) with a \\ p, there exists 
a rotation r such that the function e lp ' TX 4>(\tx — a\) = e lT p ' x <p(\x — r _1 a|) 
is permutationally symmetric. Such families describe wave packets with the 
momenta directed toward or away from the origin. 

If for d > 2 the potential V(x) is spherically symmetric, then V (0) = 
^Ay(O) • Id nxn , and therefore the eigenvalues ej of V (0) arc all equal to 
gAV(O). Thus we have 

Corollary 4.2. Let d > 2 and V satisfy Condition (VA ) and let V be spherically 
symmetric. Then L(X) restricted to permutational symmetric functions has 4 
eigenvectors or associated eigenvectors near zero: two-dimensional space with 
eigenvalue 0; and two-dimensional space with the non-zero imaginary eigenval- 
ues ±ie(X), where 

and with the eigen] 'unctions ( ), where £ and r\ are real, and permuta- 



±ir } {\) 

tion symmetric functions satisfying 



£(A) = V2j2^^o + 0(h) and n(X) = -hJ- AF(0) £ x n <& + 0(h 3 / 2 ). 



, dx n ' V d 

Besides eigenvalues, the operator L(X) may have resonances at the tips, ±iX, 
of its essential spectrum (those tips are called thresholds). Recall the notation 
a + := a if a > and = of a < 0. 

Definition 1. Let d 2. A function h is called a resonance function of L(X) 
at = ±iX if h ^ C 2 , \h(x)\ < c(x)~( d ~ 2 *>+ and h is C 2 and solves the equation 

(L(X)-n)h = 0. 
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Note that this definition implies that for d > 2 the resonance function h 
solves the equation (l+K(X))h — where K(X) is a family of compact operators 

given by K(X) := (L (A)- fx+0)- l V big (X). Here L (X) := ' ° A 1 A 



A- A 



and 



v (x)-f v h -m^f) \ 

VMa[A) -{-v h +f((^r)+2f'((^r)(^r o ) ■ 

In this paper we make the following assumptions on the point spectrum and 
resonances of the operator L(X) : 

(SA) L(X) has only 4 standard and associated eigenvectors in the permutation 
symmetric subspace. 

(SB) L(X) has no resonances at ±iX. 

The discussion and results concerning these conditions, given in |GSlj . sug- 
gested strongly that Condition (SA) is satisfied for a large class of nonlinearities 
and potentials and Condition (SB) is satisfied generically. Elsewhere we show 
this using earlier results of |CPI ICPVj . We also assume the following condition 

(FGR) Let N be the smallest positive integer such that e(X)(N + 1) > A, VA e /. 
Then ReZN+i,N < where Z m ^ n , m, n = 1, 2, • • • , are the functions of V 
and A, defined in Equation ( I14f> below (see also 

We expect that Condition (FGR) holds generically. Theorem 15.11 below 
shows that ReZ n+ i, n — if n < N. 

We expect the following is true: 
(a) if for some Ni(> N), ReZ n+ i n = for n < Ni, then ReZ Nl+ i^ Nl < 
and (b) for generic potentials/nonlinearities there exists an Ni(> N) such 
that ReZN 1+ i : N! 7^ 0. Thus Condition (FGR) could have been generalized by 
assuming that i?eZjv 1 +i,Ar 1 < for some N± > N such that ReZ n+ i, n = for 
n < N\. We took N = Ni in order not to complicate the exposition. 

The following form of RbZm+i,n 

ReZ N+hN = Im(<Ti(L(\) — (N + l)ie(A) - O)" 1 ^, F) < (9) 
for some function F depending on A and V and a± := ( ^ ^ J , is proved in 



x 1 

| BuSul IT7T1 ITY21 IT Y3I l5W4| for N = 1, and in the present paper for N = 2, 3. 
We conjecture that this formula holds for any N. [See more explicit forms of F 
for N = 1, 2, 3 in Theorem IO below]. 

Condition (FGR) is related to the Fermi Golden Rule condition which ap- 
pears whenever time-(quasi)periodic, spatially localized solutions become cou- 
pled to radiation. In the standard case it says that this coupling is effective 
in the second order (N = 1) of the perturbation theory and therefore it leads 
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to instability of such solutions. In our case these time-periodic solutions are 
stationary solutions 

'(U'"'" +,s ( ! 0'* 

of the linearized equation ^ = L(X)x and the coupling is realized through 
the nonlinearity. Since the radiation in our case is "massive"— the essential 
spectrum of L(X) has the gap (— i\,iX), A > 0, — the coupling occurs only in 
the N— th order of perturbation theory where N is given in Condition (FGR). 

The rigorous form of the Fermi Golden Rule for the linear Schrodinger equa- 
tions was introduced in |BS| . For nonlinear waves and Schrodinger equations 
the Fermi Golden Rule and the corresponding condition were introduced in |S] 
and, in the present context, in |SW31 IBuSul lB"P2l IT7T1 ITYS1 FTY5| . 

Recall that a function g S £ 2 (R d ) is permutational symmetric if 

g{x) — g(o~x) for any a € Sd 
with S n being the group of permutation of d indices and 

a(xi,x 2 ,- ■ -,x d ) := (a; ff (i),a; (T (2), ■ ■ ;x a ^). 

5 Main Theorem 

The following are two main theorems of this paper. 

Theorem 5.1. Let h be sufficiently small, the potential V be spherically sym- 
metric and the initial date ip(0) be permutational symmetric ifd>3. If f satis- 
fies the conditions (fA)-(fC) and (GWP); V satisfies the conditions (VA)-(VC) 
and (GWP); A £ I, and the spectral conditions (SA) and (SB) are satisfied, 
then there exists c, eo > such that, if 

inf{||V>o-e^M 0) £+^ 2 Vll^ 

(io) 

with < £o and z„ n = 1,2 being real, some large constant v > and 

with k = [^] + 3 if d > 3, and k = 1 if d — 1, then there exist smooth functions 
7, zi, z 2 : R+ -> K, A : M+ -> J and R : R 3 x M+ -> C sucft tftaf 

^(rs, t) = e^o Ms)ds e r 1 (t)^x + pi< j ) ^ + ip2 ^ + ( Zl + p3 )£ + l ( Z2 + Pi ) r] + R} (11) 

where p n (z,z) : C — ► M, n = 1,2,3,4, are polynomials of z := Z\ + «Z2 and z 
and of order \z\ 2 , the function R(t) can be decomposed as 

R= ]T R mn (X)z m z n + R k 

2<m+n<k<2N 

where for some large constant v {x)~ v R m _ n € £? , || (x)~ v R2NW2 J$ ™~ . 

The functions \, 7, z have the following properties 
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(A) There exists a Aoo € T such that |A(i) — Aoo| < c(l + t) ^jv as t — > oo and 
A = A m: „z m z" + Remainder (12) 



moreover 



2<m+n<2N+l 

where A TO „ = if m,n < N + I; 



(B) 

7 = r mi „z m z n + Remainder (13) 

2<m+rc<2JV+l 

where T m ^ n = if m,n < N + 1 and m ^= n; T m ^ m is real if m < N ; 
(C) \z\ < c(l + t)~2w and moreover 

z = -ie(X)z + Z m;I1 z m z n + Remainder (14) 

2<m+n<2JV+l 

2n+l,n(A) «s purely imaginary ifn<N; Z m ^ n (X) = if m,n < N + 1 and 
m =/= n + 1; moreover when N — 1,2,3 we ftave 

ReZ N+1 ^ N = Im{ai (Z(A) - (iV + l)ie(A) - O)" 1 ^, F) < (15) 

where F is some vector function such that (L(\) — i(N + l)e(\) — Q)~ 1 F = 
doRjy+i for some constant do (actually F = —doNfj+ifi with Nn+i,o 

defined in Corollary \5.4\ below), the matrix o~\ is defined as f ^ ^ ^ 
(W 

Pk (z,z) = ]T pi ™U X ) zm ~ zn ' k = 1.2,3,4. (16) 

2JV+l>TO+n>2 

Ifm, n < N then P$ n (\) and P^ n {\) are real, Pm,n.(A) and Pm, n (A) are 
purely imaginary. 

The term Remainder is bounded as 

\Remainder(t)\ < \z\ 2N+2 (t) + \\{x) -"Jfcv(t)||l + \\RN(t)\\lo , 17l 
+ \z\(t)\\(x)-^R 2N (t)\\ 2 [U) 

with v being the same in Theorem \5.1\ 

In the next we present the properties of the remainder function R in detail. 
First we define the following functions to measure the asymptotic behavior of 
the parameters and remainders: 

Z{T) :=mzx{T +t)™\z{t)\, Hi(T) := max(T + t) 1 ^ \\ P - V R N \\ nh 

JV + 1 2JV+1 

K 2 (T) :=^(T +t)Tnr||i2iv(t)||oo, ^a(T) := max(T + t)~ \\p- v R 2N (t)h 
n A {T):=^\\R N {t)\\ Hl 



(18) 
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where I := [f] + 3, T := (|z} 0) | + \z\ '\)~ x , p v = (x) v ', recall the definitions of 

zf^ and z^ and of the large constant v in Theorem 15.11 Before stating the 
main theorem we introduce the notion of admissible functions. 

Definition 2. A vector-junction u : M. d — > C 2 is admissible if the vector- 
function y p ? ^ u /ias reaZ entries. 

Theorem 5.2. The function R in Theorem 1,5.11 satisfies the equation 

Im{R, i(j) x ) = Im{R, -^-<j) X ) = Im{R, in) = Im{R, 0=0 (19) 
and £/ie function R can be decomposed as 

R= Rmn(\)z m Z n + R k (20) 

2<m+n<k<2N 

where the functions R m n ■ ^ d —> C 2 , Rk '■ U rf — * K 2 have the following proper- 
ties: 

(RA) the function R m ,n € £ 2 is admissible, and decays exponentially fast at oo 
if max{m, n} < AT; and 

i/ max{m, n} > N then the functions R m .n are of the form 

Y[(L(\)+ike(\) + 0)- n *P c <p m , n , (21) 
k 

where the function 4> m ^ n is smooth and decays exponentially fast at oo, 

< n k < N, and 2N > k > N + 1; note £/ia£ £/ie equation f \21}) makes 

k 

sense in some weighted C 2 space, see \GS£l I GSty : 

the function R k (2 < k < 2N) satisfies the equation 

f t R k = L(X)R k + M k (z,z)R k + N min{k ^ N} (R min{k ^ N} ,z,z) + F k {z,z), 

(22) 

where 

(1) F k (z,z) = 0(|z| fe+1 ) is a polynomial of z and z with X-function-valued 
coefficients, and each coefficient can be written as the sum of functions of 
the form /W\). 

(2) M k (z,z) is an operator defined by 

M k {z,z) :=jP c J + \P cX +A k (z,z), 
where A k (z, z) is a 2 x 2 matrix bounded in the matrix norm as 
\A k (z,z)\<c\z\ 2 e-^. 
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(3) Nn{Rn, z, z) satisfies the estimates 

\\n n (r n ,z,z)\\i + \\Nn(Rn,z,z)\\ H 2 < c(T +t)-^[z 2 nln 4 2 +n 5 2 Til], 

(23) 



||(-A + l) 2 iV w (i? w ,z,z)|| 1 + ||(-A + l) 2 7V j v(i? A ,,z ! z)|| 2 

< c{T + t)-^{z 2 nlnl + z 2 nii^nl + i$ni). 1 ' 

Proof of Theorems |Oand 15721 

By plugging Equation ( lllfl into Equation ( 2J we have 

A<^ + d t \pi<f>$] + ii£ + \[zi + P3 ]Ca 
dt\p2<l> X } + z 2 r) + \(z 2 + Pi)-q\ + (zi +p 3 )e(X)i] 

(25) 

with a=:(l), Ri := JteH, R 2 := JiV(i?, P , ,) := ( 

and 

ImN(R,p, z) := /(|0 A + A + il 2 \ 2 )l 2 - f{{<j> x f)h 
ReN(R,p, z) := [f(\<? + h + il 2 \ 2 ) - f {{<t> x f)]{4> X + h) - 2/'((0 A ) 2 )(0 A ) 2 / 1 
with 

h :=Pi^a + (*i +Pa)Z + Ri, h=P2<\> X + {z 2 +p±)t) + R 2 . 

From Equations ( I25JI and the orthogonality condition ( I19JI we obtain equa- 
tions for Z\, z 2 , A, 7 

[j t (z 1 +p 3 )-e(\)(z 2 +p 4 m, V )-(ImN(R,p,z), V )=F 1 ; (26) 
[j t {z 2 +Pi) +e(X)( Zl + P3 M,v) + (ReN(R,p,z),0 = F 2 ; (27) 



7 + d tP2 - Pl + 1 (ReN(R, p, z) , cp x ) = F 3 ; (28) 



^ + dtPi- 7 -^— I r(ImN(R,p,z), ( j ) x }^F il (29) 



where the scalar functions F n , n = 1, 2, 3, 4, arc defined as 

Fi := j((z 2 +P4)r),v) - x (i z i + P3)Z\,v) - 7(^2,?/) + \{Ri,7}\); 
F 2 := -j((zi +p 3 %0 - \((z 2 +P*)rix,0+j(RuO + A(i? 2 ,6>; 



1 
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Moreover using the fact that P c ^i? = j^R — XP C \R we rewrite Equation ( I25fl 

as 

j t R-L(X)R-P c JN(R,p,z) = g (30) 

where the function 



Q := -XPcxR + jP c JR + P c 



7(f2 +Pi)V- ^Pidi(j) - AOi + p 3 )£,x 
-Xp 2 (j)\ _ A(z 2 +Pi)rj X 



In the next we sketch the proof of the theorem except Equation ( I15|l . We 
start with proving the expansion of R and the expressions for z, A, 7 and 
Pi, I = 1,2,3,4, by induction on k = to + n, k = 2, 3, • ■ •, with the main tool 
Lemma 15.31 below. For the detail of finding R m ^ n , A m „, T rn ni Z m n we refer 

to |(tS2I IHS3| . For k — 2, by Lemma 15.31 we have that iNm]n, defined in 
( I32|) below, is admissible (Note Nm] n , m + n = 2, does not depend on any 
Rm',n' , m + n > 2 by Lemma 15.31 and the term C? has no contribution for terms 
of order \z\ 2 ). Thus R m ^ n = ~[L(X) + i(m - n)e(X)}- 1 N$ n , m + n = 2, is 
admissible if N > 1. Moreover by the properties of Z m . ni A m „ and r mjra in the 
equations for A, 7, i = ii + 12:2 in (EJ-(E}> the equations ( l2"o| -fl2"§ )l and the 
observation that the term Fj, j — 1, 2, 3, 4 has no terms of order \z\ 2 we have if 
TO + n = 2 and m 7^ n + 1 then 

H(m - n)e(A)P$, - e(A)P^l] + i[-i(m - n)e(X)P^ n + e(X)P$ nl 
[(ImN$ n , r,) ~ i(ReN$n, 0] 

and if to' + n' = 2 and to' ^ n.' then 



-i(rnf - n')e(\)P<$ in , - P«„, = --—(ReN^ 

{9^ Ax) 

-i(m> - n')e(A)P™ = -^(I m N™<p) 



where ( ^ ™'(2) ) : = Nm]n- We observe the solutions satisfying the condi- 



-iZeiV, 



tion Pm,n = Pn,m, = 1, 2, 3, 4, to the equations above exist, and by the prop- 
erties of Nm,n we obtain Pm,n(A) and Pm,n(A) are real, Pn,n(A) and Pm,li(A) 
are purely imaginary. Thus we finish the proof of the step of to + n = 2. For 
the step m + n = 3, P^j; R m > tn >, ml + n' = 2 and Lemma 15.31 enable us to 

determine Pm]m R m ,n, m + n = 3 and prove all the properties in the theorem. 
By using the procedures above repeatedly we finish the proof. 
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By using the same techniques in |GS21 IGS3j and Equation ( l2*2*)l we prove 
the decay of Rn and R2N in appropriate norms; and by f I12|l there exists a new 
parameter [3 — z + 0(\z\ 2 ) such that 

$ = -ie(X)[3 + E Y n (3 n+1 fi n + Remainder, (31) 

l<n<AT 

where Y n is purely imaginary if n < N and ReY^ = ReZ^+i,Ni consequently 
we have 

- — \z\ 2 = ReZ N+hN \z\ 2{N+1) + \z\Remainder 

which implies the decay of f3, hence y, by the assumption ReZjsr-\-i t N < and 
the estimates on \z\Remainder; and similar as in |OS3| by the expressions for 
A in ( I12JI there exist c m)ra (A) such that 

— [A— E] c rnn [X)z m z n \= Remainder 
2N+X>m+n>2 

with the the term Remainder satisfies ( 1170 , consequently A is convergent due 
to the decay of z and the fact that Remainder(t) is integrable at 00. 
The proof of Equation (EJ is m Sections [7|and below. 

□ 

In the proof we use the following lemma. 

Lemma 5.3. (A) If Gi is a vector-function from M. d to C 2 such that iG\ is 
admissible and G\ £ C 2 , then the vector function (L(X) + ifi)~ 1 P c Gi is 
admissible for any \i G (—A, A). 

(B) Let Rk, R m ,n, Pv be the same as those in ( \2U\) and respectively and 
Pv= E <^?n(A)^ n , v = 1,2,3,4 

m+n>2 

with gm,n(A) and qin,n(X) are real, qm]n(X) and qm,n(X) are purely imagi- 
nary if m, n< min{fc, N}. Then JN(R,p, z) can be expanded as 

JN{R,p,z) = E z m z n N^l(X)+A k (z,z)R k 

2<m+n<2N 

(32) 

+^min{fe,AT}(-Rmin{fe,JV},P, z) + Remainderl 

where iNm}n(X) : R 3 — > C 2 is admissible if m,n < N and only depends on 

Rm' .ri Qrn 1 n' w ^ n m' < m,n' < n and \m' — m\ + \n' — n\ 7^ 0; the 2x2- 
matrix function Ak(z,z) is bounded in the matrix norm as \Ak(z,z)\ < 
c|z||e _eo ' a: '; Nj(Rj,p, z), j < N, contains all the nonlinear terms in Rj , 
Nn(Rn ,p, z) admits the same estimates as in f \2<S\) ( \24\) ; and the term 
Remainderi has the estimate 

\Remainder x \ < c\z\ 2N+1 e- e °^ . (33) 
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Proof. The first part is copied from |HS2llHFvp . 

For the second part we prove the expansion for JN(R,p = 0, z) in [GS21 
ICjS3| . The proof is almost identical to that, thus omitted. □ 

The following corollary will be used later, whose proof is obvious by f 
thus omitted. 

Corollary 5.4. Let k = 2N in then 

2N 

JN(R,p,z)= ^2 N m>n (X)z m z n + N N (R N ,p,z) + Remainderl 

m-\-n=2 

with iN min admissible for m,n < N; Nn(Rn,p, z) and Remainderl have the 
same estimates as above. 

6 The Effective Equations for i, A, 7 and R 

The equations for the parameters i, A, 7 and the function R have some terms 
having no role on determining RcZn+i,n- In this subsection we separate the 
effective parts from the "useless" parts. We will use the following lemma to 
show that the terms on the right hand side of Equations ( I26|) -f Tlty and ( I30|) 
are useless. 

Lemma 6.1. If we expand the scalar functions F1+1F2, F3, F4 and the function 
Q in z and z, then we have 

F 1 +iF 2 = K^JX)z m z n + Remainder (34) 

2<TO+n<2JV+l 

F 3 = K ™,nW zm z n + Remainder 

2<m+n<2N 

F 4 = Y K ™,nW zm z n + Remainder 

2<m+n<2N 

G = Y Gm,n( X )z mz ~ n + GN{R N ,z, z) + Remainderl (35) 

2<m+n<2N 

where if m,n < N + 1, then the coefficients Km}n, Km^n are purely imagi- 
nary, Km]n is real, the function iQ m ,n is admissible; Km,n = and Q rn ,n = 
if (m, n) = (0, 1) or (I, 0) with I < 2N; moveover the function Remainderl sat- 
isfies the estimate ( \S3\) . the term GN(Rn,z,z) admits the same estimates as 
N N (R N ,z,z) in 

Proof. By the expansions for the scalar functions A, 7 and the function R in 
Theorems 15.11 and 15.21 we prove the lemma by direct computations. □ 
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7 Proof of Equation ( [T5ft when N = 2 

In this section we prove Equation f ll5fl when N = 2. For consideration of space 
we only prove the case of the nonlinearity f(x) = x. One can show that the 
results hold for more general nonlincarities, but the computations are much 
messier. The main result is 



Theorem 7.1. If N = 2 and the nonlinearity f(x) — x in (^\), th 



en 



ReZ 3 . 2 = -^—Im{ ai {L{X) - Zie{\) - 0)- l P c N 3 , o ,N 3 , ), (36) 

where, recall the definition of N 3> q from Corollary \5.4\ 
The proof is in subsection 17.21 

7.1 the coefficient of z 3 z 2 in the expansion of (ImN(R,p, z),r)}- 
t(ReN(R,p,z),0 

Suppose that the expansion of (ImN(R,p, z),rf) — i{ReN(R,p, z),£), the part 
on the left hand side of Equations ( I26[) and ( 1271) , in z and z is 

(ImN(R,p,z),Tj) — i(ReN(R,p,z),£) = X m . n z m z n + Remainder. 

2<m+n<5 

In the next we compute X 3y 2- Before starting the proof we define the following 
functions and constants: 

K ._( -P$(<P X n + iP£U X <P X xv\ K ._ ( i<P X V, \ R 
K ._ 1 ( + ifr, \ K . ( -^[P 2 ( J - %Pj$] \ 

Ks - 8 +e J'^-l -i*vp« + wep^ ) 



and 



D 1 (X) : =i[P 3 ( > \ ) (^(r ? 2 -3a,^)-iP3 ( ,i ) <2^ A ,^)], 
D 2 (A) := -2i{axR 3fl , #i + K 2 + 3K 3 + K 4 ), 



D 3 (X) := (P 3> i, 



-i((f> x ri 2 - 3(/> A £ 2 ) 
20 A £r? 



and 



Z? 4 := pJJ [-3i<^ 2 , + iifcr), yy)] + P 8 <g (20\ 2 , 



and the matrix a\ := I _ ^ ) , recall the definitions of Pm,n, I — 1,2,3,4 
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in Theorems 15.11 The main result is 
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Proposition 7.2. 

4 4 



n— 1 n— 1 



The proof is at the end of this section. In the proof the following observation 
will be used 



Lemma 7.3. Pm,n = Prt.m for any m, n and k; If k = 1,3 and m, n < 2 then 
Pm}n is real; if k — 2, 4 and m,n < 2 then Pm^n is purely imaginary. The 
constants Pm,n have the follow expressions 

^°-4 £ (A)(^,^) (0 ^ 0) ' z2 £ (A) + 8ie(A)<*\#) ' 



o, p(3) pW _ 1 ip ( -icf>W . 
" 2ZP3 - 1 _ ^ " <AK^> ( 3 '°' I <^7 1 + ' 



jl ~~ 2e(A)(5'(A) \ — «(</> A ) 2 £ ) )+U 

,x ~~ _ i2e(A) ~ 2e(A)<5'(A) ' 3 ' ' y P + %U ' 

U stands for real ( different) constants throughout the rest of the paper. Ifm,n< 
3 and \m — 3| + \n — 3| ^ t/ien 

iZm,» = -(L(X) + i{m - n)e(X) + O)- 1 ^™,™ - Q m ,n\- (37) 



(k) (k) 

Proof. The property P m ' n — Pn,m following from the fact that A, 7, z\, z 2 are 

(k) ' 

real functions. The fact Pm'n being real or purely imaginary is based on the 
properties of N m ^ n in Corollary 15.41 By Equations f!26 |) -f[2l? jl and the properties 
of A, 7, i in Theorem 15.11 we compute to get the constants P m .n- 

The formula of R m .n follows from the expansion, we refer to our paper 

321 |QS3|. □ 



We have the important remark for —0 symbol in the expression of R n 
Remark 4. If \m — n| < 2 then 

(L(A) + i(m - n)e(X) + = (L(A) + i{m - n)e(A))- 1 . 
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Proof of Proposition 17. 21 For the fact ReX^ ^ = Yln=i R e D n {X) we only 
use Lemma 17.31 the fact that i?„ i; „ is admissible if m, n < N = 2 in Theorem 
15.21 and direct computation, thus omit the detail. 

For the fact on the right hand side we transform the expression for D n , n = 
1,2,3,4. 

First we prove the D\ part. By the expressions of P 2 ,o, ^3 n ^20i ^3 l 
and the properties of Pm,n m Lemma 17.31 we have 

y-3£ a ),$> + 



ReD 1 = Rei[P$ {4> x {r? - 3£ 2 ), + 2P 3 (2 1 ) (-i^, A ) 



= -Re{R 3t0 , AioxKx). 

For D3 we compute P^.o and ^3,1 first. Recall the forms of R m ,n in Lemma 
17.31 and the properties of Gm.n in ( I35fl . First we have 

P 2 ,o = ~[L(X) + 2^e(A)]- 1 P c iV 2 , (38) 
3^A^2 _|_ 0X^2 J ■ The function P 3i i has two parts 

i? 3 ,i = ~[L(X) + 2^e(A)]- 1 P c (G 1 + K 5 ) (39) 

where the vector function iG\ is admissible, hence — [L(X) + 2ie(A)]~ 1 P c Gi is 
admissible (thus " useless" ) by Lemma 15.31 and the vector function K§ is 

which is part of iV^i. By the observations that ai(L(X))* rj\ = L(X), a\ = — 1 
and the definition of K2 we have 

PeD 3 = -Pe(P 3 , ,4iaiX 2 ). 

D4 admits the form ReD^X) = 4Im(<JiR3fi, K4), which follows from some 
manipulation on the expressions of P^q and Pj™"*. The proof is tedious, but 
not hard, thus omitted. 

Collecting the computation above and recalling the form of D2 we have 



J2 ReD n {X) = 6/m( < 7 1 P 3) o, J2 K ^ 



n—l n—1 

Thus the proof is complete. 

□ 
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7.2 Proof of Theorem [7TT1 

Proof. Recall that X 3 2 is the coefficient of z i z 1 in the expansion of the term 
(ImN(R,p,z),ri) -i(ReN(R,p,z),£). First by Equations (H>J), (EJ and 
we have ReZ 3 2 — 7r^y-R e ^3,2, then by Proposition 17.21 

4 4 

i?eX 3 , 2 = ^i?e£„(A) = 6Jm<(Ti723,o,53-Kn>, (41) 

n— 1 n— 1 

where, recall the definitions of if n before ( 1361) . 

By (ED and the properties of Q in (ESI for TV = 2 

#3,0 - -(i(A) + 3ie(A) + Q) _1 P c iV3 i o 
and we compute to get 



7V 3 ,o = - £ ^ 



Consequently we have 



ReZ 3>2 = ^J^ReD^X) 

n=l 

= ^Im<£n(L(A) + 3ze(A) + Oj-^cJVa.o, A^o) 
= ^J m ( CTl (i(A) - 3ze(A) - 0)- 1 P c JV 3l o,iV3,o>. 

Thus the proof is complete. □ 

8 Proof of Equation ( CGD) when N = 3 

The proof of N = 3 is more involved, but the ideas are the same. Thus in this 
section we just state the key steps and only prove the case f(x) = x and the 
potential V is spherical symmetric. The following is the main result. Recall 

that G\ = ( j 0^1 an< ^ ^ 6 a * e ^ m ^^ on °^ ^ 4 -° f rom Corollary 15.41 

Theorem 8.1. If N = 3 and the nonlinearity f(x) = x and the potential V is 
spherical symmetric in f \J\) . then the constant Z^ 3 satisfies the equation 

ReZ±. 3 := -^-I m ( ai (L(X) - 4ie(\) - O)" 1 ^, N ifi ) 

where, recall the definition of N^.q from Corollary \5.4\ (for N=3). 



In the next wc state the following observations which will simplify our com- 
putation. We defined two spaces 

$ = {g: M. d -► C| g(-x) = g(x)}; AS = {g : R d -» C| ff (-a?) = -g(x)}. 

Obviously S ± AS and <jl A ,^ £ 5 and £,77 e AS (see Corollary |OJ. Recall 
that N m ^ n and 5m,n are defined in Corollary I5.4l and Equation ( 13 5(1 respectively. 
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Lemma 8.2. If the potential V is spherical symmetric, then R m ,m A'm.n € § 
if m + n= 2, 4, 6; R m ,n, JV mi „ € AS if m + n = 3,5. Moreover 

p(3) = p(4) = p(l) = p(2) = n / 4 ox 
m,n m,n m 1 ,n' m' ,n' \*"J 

if m + n is even or ml + n' is odd. If m,n < 3 i/ien Pm.n, -Pm,n are real and 
Pm^n, Pm,n are purely imaginary, and 

P$n = Pk% ifk=l, 3, P«„ = -P(% if I = 2, 4. (43) 

If m,n < 3 and \m — 3| + \n — 3| 7^ i/ien 

iim.n = -(L(\) + i(m - n)e(X) + O) -1 ^,,, - Q m , n \. 
(k) 

Proof. The facts Pm'n being real or purely imaginary follow from similar dis- 
cussion in Lemma 17.31 The rest follows from direct computation based on the 
properties of N m<n in Corollary 15.41 and the observations made on £, 77, 4> x , <j)\ 
before the lemma, thus we omit the detail. □ 

8.1 The Explicit Forms of i? 4 and A^ j0 

We denote the coefficients of z m z n in the expansion of ImN(R,p, z), (ReN(R,p, z)) 
as ImN m ^ n , (ReN m n ). The following lemma is the main result of this section. 

Lemma 8.3. 

R 4fi = -(L(\)+4ie(\) + 0)- 1 N 4 . (44) 
with Ni,o given explicitly below. 

The proof follows from the definition of R^o in Lemma 18.21 the property 
£4,0 = in (ESJ) for N = 3. 

In the next we compute N 40 = ( ^t, tf ) explicitly. As usual we denote 

^ — RelM^o J 

( Rnin \ (k) 
R m n — I ^ n J and recall that for some m, n, k, Pm.n = from Lemma 18.21 

We compute to get 

3 3 
ImN 4fi = HP 1 ^eiV 4 , = H n } ( 45 ) 

n— 1 n — 1 

(k) 

with the functions Tin , n = 1, 2, 3, ft = 1, 2 3 defined as 

+2^ x p^r% + \e^l - 3 4 v 2 r% 
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+ 2(0 A ) 2 p 2 > 2 (2 o ) - 

Hf := 30 A £P« - ^ A r?Pg. 

8.2 The Explicit Form of Z 4 3 in ( ED 

The following lemma is the main result of this section. 
Lemma 8.4. 

ReZi. 3 = --^—Re[E 1 +E 2 + E 3 ] (46) 
where the constants E n , n = 1,2,3, are defined below. 

Proof. By Lemma 16.11 we have that the right hand sides of ( I26|) and ( I27|) are 

useless in determining the real part of Z4.3. Thus ReZ^ j, = Re\{ImN^, n) — 

i{ReNi t 3, £}] by the fact z = z\ + iz2- What is left is just direct computation 

(k) 

while deserting useless terms by using the property of Pm, n , k — 1,2,3,4 in 
Lemma 18.21 and using the fact R m ,m rn,n < 3 are admissible. □ 

In the next we give the definitions of E n , n = 1,2, 3, and the various related 
terms. 

Ei := (#4,0, ( J )> + (#4,2, ( p e ))-4i(R 4 ,i, ( ^ ))-P$G 5 -P$G e ; 

E2 '■= P^\{Gl + G5) + P4 i(G2 + Gq) + P42G3 + #4,2 G4,; 

E 3 := (i?4,J \)+4i(R i<u ( ^ 

T x := -20 A ry 2 P 3 (4 o ) - iS^Pg - 2^J?g> - *6^i$o + ^V P 2,o 
-20 A ^P 2 (2 O ) - 2^Pg + " S^Pgo - 3i^ 2 P$; 

T 2 := -2^P 3 ( g - 2i^r,P^ - 2<^P$ - 2# A £Pjg " 2<^P 2 (1 ) 
-2^P« + 3^ 2 i45 - i?R<» + 3i*VP$ - i^Pg; 

F 3 := -2(0 A ) 2 r ? P 2 (2 o ) - 20V? 2 2 6' + f ?V 

JP 4 := -20V^ 2 ( ,o - 2 ^ Aj ?4!o " zd 2 V 
-2iPg(^) 2 C+| ?7 3 -2#^45; 
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T 5 := #V - z30 A £ 2 ; 
:= -2^fo 

Gi := -<2(^) 2 £?P$) - 2(<^Pg> - f (£V> 
+<6i^0fcVo ( ,a ) + I <e 4 ) + 6i(<^ 2 pg>; 

G 2 := 2iP ( $((^) 2 £r 7 ) - (20V^ 2 )P O ( ,2 ) - 2faVi$) 

+ §<^)-§<^)+i(2^J@; 

G 3 :=-i<^7j 2 )+3^ A ^ 2 ); 
G 4 := -((«^) 2 £t?); 
G 5 := 12z(^ 2 ,P«)-4(^,P$); 
G 6 := -4(0 A r, 2 , i?^) + P$) 

where we use (gi), (.91,32) to denote J gi(x)dx, J g\g 2 dx respectively for any 
function g\ and g 2 . 

8.3 Expansion of E n , n= 1,2,3 

(k) 

In this section we first relate the terms i?4,i, i?4,2, P\ \i k = 1,2,3,4, to the 
function R^q so that ReZ^^ = Re(R4.o, 1 $) for some function "P. Then after 
some manipulation we relate the function iff to A^o, hence to i?4 i o- Recall that 
the relation between P40 and A^o in Lemma 18.31 

The main result is the following proposition whose proof is in the later part 
of this section. 

Proposition 8.5. 

Re[E 2 ] = Pe[P 4 ,o + £4,1 + P^M + PftN] , (47) 

and 

Re[Pj$Yi + P S Y ^ = Pe[(P$, -QiH^) + (R®, -6iH^)}, (48) 

Re[Ei + P 4 ,o] = ReiR^SiH^ -2iH ( 2 ) -2iH ( 3 } ) (m 
+Re(R^l-8in[ 1] -2iH ( 2 1) -2iH ( 3 1) ), [ ' 

Re[E 3 + £4,1] = Pe[(p£>, ~6iH { 3 2) ) + (P$, -WH^)] (50) 
where the functions Y\, Y 2 , E^ and E41 are defined below. 
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In the decomposition of E<i above we used the following constants 
P 4>0 := -MRU tfflpg + Mf )pW) - 4*<pg, M«P,g + Mf >P «) 

P M := -4i(pg, M 3 (1) P ( 5 + Mi 2 )p«) - 4i(J®, Mi 1 )^ + M^pg), 
F a := -Wf'Pj'o 5 ~ W^P-ft + d + G 5 , 



and 



Y 2 := -4iW, (1) P (2) - 4ziy 9 (2) P 9 (1 n } + G 2 + G 6 



Recall the notation 5 (A) := j^m^) ^ A )- Before starting proving the propo- 
sition, we compute the explicit form of P^l , P4 2 : 

-2ie(A)(5'(A)P 4 (1 2 ) 

= (fl«^f») + (P$,M«) + (Pi:i,M«) + (Pi 2 l,MW) (51) 
+P^W^+P^W^+iU. 

where, {/ is a real constant and ill includes all the contributions from F4 (see 
Lemma 16. 111 . 

M 1} :=2(^) 3 P 2 ( S+2(^) 2 PS-5^ 
M« := 2(0 A ) 2 ^P 2 ( , 1 O ) + 2(^) 2 i?« + \^ - \<t>W 
M« := -i(^f V := (0 A ) 2 £, 

Wf> :=z((^) 2 C7?), M/ 2 (1) :=((0 A ) 2 ^), 



~2ie(A)(5'(A)P 4 (2 2 ) = (^f + S)!^^ 



-(<^M 2) + ^) + «l,M 2) + ^S) (52) 



where J7 is a real constant including all the contributions from F3 (see Lemma 

EH, 

A*j» := -6^) 2 P« - 6^*$ - f ^ + I^^ 2 . 

X 2 2) := -2(0 A ) 2 A P 2 (2 O } - 20 A A P 2 2 6 ' + i^fr, 

M| 2) := -30 A A £, M { 4 2) := z0 A A 77, 
In the following lemma we relate Y"x and Y2 to P3 q and p{ 4 J . 
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Lemma 8.6. 

Yi = 6ie(A) (£, rj) [3»P 3 ( J - P 3 (3 } ] , F 2 = 6ze(A) (£, ry) [3iP 3 (3) + pg]. (53) 
Proof. Using Equations ( ( 113 to obtain 

-8e(A)<£, ^Pg = -3i(/miV 3)0 , 77) - (PeA^, £) (54) 

and 

8e(A)(e,ry)P 3 (4 ) = -3i{ReN 3<0 ,0 + (ImN 3>0 , V ). (55) 
On the other hand we observe 

Yi = |i<e 4 ) " f <£V> + 6((0 A ) 2 ^)P 2 ( .o ) + 6<<^?Pg!> + 18 i (^^ 2 P^) 
+ 18 i (0V^ 2 )P 2 (1 o ) 
= 6i<JteJV 3>0 ,£) 

and 

F 2 = |(?7 4 > - Kerf) 6iP a ( J<(^) 2 ^) - 6<0V^3) " 6i<^J2$> 
= -6i(ImN 3 .o,ri). 

These observations together with ( I54|l and ( I55|l imply ( I53fl . □ 

Proof of Equations f l47[l and (I48[) By the definition of P 2 it is sufficient 
to prove ( I47f) by proving 

Re [p[ 1] 2 G 3 + P$ G 4 ] = Re[E 4t0 + E 4A +P${Y 1 -G 1 -G 5 ) + P 4 (4) (Y 2 -G 2 -G e )}. 

(56) 

First we observe that 

Gs " 2^A) G4 = W p *% G 4 = -8e(A)5'(A)P 2 W. 

Lemma 18.21 implies that P^q and P 2 q are real and purely imaginary respec- 
tively, this fact together with Equations ( I51fl ( I52H implies ( I56|) . The compu- 
tation is straightforward, but tedious, hence we omit the detail. 

Now we turn to ( I48|) . We start with analyzing P| 3 { and P 4 4 j . Equations 
( 1261) ( I27|l and the properties of the equation for z in ( I14fl imply 

-8e(A)(£,7 7 )P 4 (3) = SiilmNi^r]) - (ReN^i, £) + U (57) 

and 

8e(A)<£, 7 ? >P 4 (4 1 ) = -3t<iteJV4,i,0 + {ImN^rf) + iU (58) 

where f stands for real (different) constants including all the contribution from 
the terms Pi and P 2 (see Lemma 16. If) . Moreover we compute to get 

( ImN^ \ ( i^vBgl + fongl \ 
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where U± is a vector function such that iU± is admissible. Thus we have 
(ReN itl ,$ = (i? 4 ,o, ( Wf/ 

and 

(ImN^r,) = (P 4 ,o, ( 'if^ ^)+iU 

Collecting the computation above and using the properties of Y\ and Y2 in 
f 153)1 and the facts that P$ = P$, P$ = -P$ in Lemma |Owe have (EEJ). 

□ 

Proof of Equation ( 149(1 We start with rewriting (i?4,2, f ^ 5 )) to relate 



6 



the function ( J to R2 q. By the definition of R4 2 in Lemma 18.21 we have 

V J-6 J 

Ra.2 = -(i(A) + 2ze(A))- 1 P c [iV4,2 + E/i] 

where C/i is a vector function such that iU\ is admissible and includes all the 
contributions of Q in (see Lemma IB.ll for N = 3). 

Recall that 01 = ( ^ 1 ]. The facts <7i(L(A))*<7i = L(A) and a\ = -1 



imply 



(i? 4 ,2, ^ J) = <iV 4 ,2,-cri(L(A) + 2ie(A)) _1 F c CTi ^ ^ J)+it? 

where {/ is a real constant generated by JJ\ above (see Lemma 15. 311 . Moreover 
we observe that 

#2,0 = -(L(\) + 2ie(\))- 1 P c N 2fl 

= i(L(X)+2ie(X))- l P c a 1 (^ 5 



consequently 



(R4.2, ( J 5 )> = -4i(iY 4: 2,aii?2,o) +iU. 



In the next we decompose N^2 into various terms. By direct computation 
we have 

icj) x r) </> A £ 



N 4 .2 - MR 4fi + ^ _ j i? 4 ,i 

z^r? ^ \ ( Pi? V „ 

where C/i is a vector function with iU\ being admissible, M. is an 2 x 2 matrix 
defined as 

/ 2<^P<g + |^ + 2(^) 2 P$, ~ I'? 2 + 2^^P ( j 1 2 ) + 2<^<2 \ 

V + S^^Pq^ - + 6^4% 2<^pg + 2(^) 2 P (2 2 ) + T 
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The definitions of E\ after ( yields 

Ei = <i? 4 ,o, ( J ]) -4i(A4E 4) o, ^1^2,0) +iU 



which together with the definition of i?4,o, Lemma I8.2l and the fact R m ,n, m,n < 
3 is admissible implies Equation ( I49f) . The computation is straightforward, 
hence we omit the detail. 



□ 



Proof of Equation ( I50[) The definitions of E3 and E^,\ imply 

E 3 +£4,1 = {IU,i,K)+iU 

with [/ being a real constant and 

K ( F 3 -4in 1 + 4i[MPP$+MWp$] \ 
4z0 2 + 4^[Mi 1) P (2 2 ) + P$] J ' 

In the next we relate the function K to i?3,o- 
First we have 

i? 4 ,i = -(L(X) + 3ie(A))- 1 F c [^ 4 ,i + Ui] 

where as usual U\ is a vector function such that iU\ is admissible and includes 
all the contributions from Q\ from ( I3U|) . On the other hand by Lemma 18 .21 and 
direct computation we have — 6^01^3,0 = K, thus 

£3 + £4,1 = 6»<i24,l,(TiJV3,o) 

= 6i{Ni i i,—(L(X)* + 3 l e(A))- 1 P c V 1 7V 3 , > + iU 

with U being a real constant generated by U% above (see Lemma I5.3|l . By the 
facts that ai{L(\))*a\ = L(X), o\ = — 1, — aiP*a% = P c and the definition 
R&,o = -(L(X) + 3ie(A))- 1 F c iV3,o we have 

E 3 + E 41 = -6i{N AA ,aiR 3t0 ). 

By the explicit formula for iV^i in ( I59|l . we obtain 



(2) 
3 



£J 3 + £41 = <i?4,o, - 6t I ^fi) J)+if7 



for some real constant {/. The proof is complete. 



□ 
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8.4 Proof of Theorem HO 

Proof. Recall the explicit form of Z43 and R4 in f I46|) . By Equations ( I47|l - 
( EDI and the observation that N 4 , = ( H n 5 > _ XI W » } J in ( E3> we 

\ n=l n=l / 

have 

ReZ 4t3 = J ^Im(a 1 (L(X)+Ate(X) + 0)- 1 N 4fi ,N 4fi ) 
= T ^Im(a 1 (L(X)-Ate(X)-0)- 1 N 4fi ,N 4fi } 

to complete the proof. □ 
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